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ON THE FORMATION OF BROWN DWARFS 

Ing-Guey Jiang 1 , G. Laughlin 2 and D.N.C. Lin 2 
RESUMEN 

Las propiedades de las enanas cafes presentan problemas para la teorfa de la formation estelar. Como sus masas 
son mucho menores que la masa de Jeans de las nubes interestelares, las enanas cafes se forman probablemente 
por fragmentation secundaria, y no por cl colapso directo del nuclco de la nubc molecular. Para evitar la 
acrecion de masa posterior, las jovenes enanas cafes deben salir de las regiones de alta densidad donde se 
formaron. Proponemos que las enanas cafes se forman en las regiones externas, opticamente delgadas, de los 
discos circunbinarios. Las interacciones dinamicas subsecuentes con sus binarias causan que las enanas cafes 
sean dispersadas a grandes distancias, o que escapen del sistema a baja velocidad. 

ABSTRACT 

The observational properties of brown dwarfs pose challenges to the theory of star formation. Because their 
mass is much smaller than the typical Jeans mass of interstellar clouds, brown dwarfs are most likely formed 
through secondary fragmentation processes, rather than through the direct collapse of a molecular cloud core. 
In order to prevent substantial post-formation mass accretion, young brown dwarfs must leave the high density 
formation regions in which they form. We propose here that brown dwarfs are formed in the optically thin 
outer regions of circumbinary disks. Through post-formation dynamical interaction with their host binary 
stars, young brown dwarfs are either scattered to large distance or removed, with modest speed, from their 
cradles. 
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1. INTRODUCTION 

Brown dwarfs are entities with mass below that 
require for hydrogen burning to ignite (< 0.075Mq) 
and above that associated with gaseous giant planets 
(~ a few Jupiter mass, Mj). Although the existence 
of brown dwarfs was proposed by Kumar in 1963, 
their cool dim nature consigned them to a strictly 
theoretical status for more than three decades. Re- 
cently, however, improved observational capabilities 
have led to the discovery of many brown dwarfs, 
prompting a renaissance in our understanding of 
these objects. Because the mass of brown dwarfs 
is much smaller than the usual Jeans mass for a typ- 
ical molecular cloud, Lin et al. (1998) claimed that 
the encounter between two protostellar discs might 
increase the Jeans mass locally. Reipurth & Clarke 
(2001) suggested that brown dwarfs are substellar 
objects because they have been ejected from new- 
born multiple systems. They use a simple model of 
timescales to show that this could happen. Bate et 
al. (2002) used a smoothed particle hydrodynam- 
ics code to show that brown dwarfs could be formed 
through the collapse and fragmentation of a turbu- 
lent molecular cloud and thus confirmed what was 
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suggested by Reipurth & Clarke (2001). 

However, if all brown dwarfs have to be ejected 
to avoid accreting too much mass, one might ask 
how it would be possible that they are wide com- 
panions of solar-type stars sometimes. In this paper, 
we propose a formation scenario for brown dwarfs 
which provides a natural explanation for the cur- 
rent observational situation. We suggest that brown 
dwarfs could be formed through disc fragmentation 
and we study the "escape zone" where the brown 
dwarfs could be ejected and become a field star. 

In principle, we demonstrate that in some cases, 
the brown dwarfs get ejected but in other cases, they 
could become long-period companion. We also study 
the criteria that the brown dwarf pair can survive. 
However, we do not intend to study the problem of 
brown dwarf desert here. 

2. STABILITY OF BROWN DWARF 
COMPANION AROUND BINARY STARS 

We now investigate the orbital stability and evo- 
lution of brown dwarfs formed in circumbinary rings. 
These companions are perturbed by the tidal dis- 
turbance of the binary star's gravitational poten- 
tial. For computational convenience, we assume 
these newly fragmented brown dwarfs quickly be- 
come centrally condensed and the residual gas does 
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not contribute significantly to the gravitational po- 
tential such that the dynamics of the system may be 
described by a few-body approximation. 

In this section, we further simplify the interaction 
procedure to a three-body (the host binary star plus 
a brown dwarf) problem. The mass of each brown- 
dwarf fragment is sufficiently small that they do not 
significantly perturb each other on the short term of 
several thousand binary periods. In order to make 
a direct comparison with some existing results, we 
first treat the brown dwarfs as massless particles. 
But, in general, we carry out full 3-body integration 
in which contribution due to the finite mass of the 
brown dwarf is included. 

We consider a range of ratio (fj, = Mi/(Mi+M 2 )) 
of masses (Mi and M2) for the two components of 
the binary stars. Following the approach by Hol- 
man & Wiegert (1999), we consider, for the host 
binary system, a range of orbital eccentricity (e»). 
The semi-major axis of the binary is set to be unity 
such that all other length scales are scaled with its 
physical value. We also adopt G(M\ + M2) = 1 such 
that the binary orbital period is 2*7r. 

Since we assume they are formed in circumbi- 
nary rings, we consider brown dwarfs with orbital 
semi major axis larger than that of the binary sys- 
tem. At the onset of the computation, all three stars 
are located at their apocenter with respect to their 
common center of mass. It is possible that the frag- 
ments may have a range of orbital eccentricity (e^). 
Here, we consider two limiting eccentricity for the 
brown dwarf (ef, = and 0.4). We also set the ratio 
(p,b) of brown dwarf's mass to that of the binary to 
be 0.05 and we assume brown dwarfs rotate in the 
same direction as the orbit of the binary secondary. 

2.1. Ejection criteria 

We choose a range of initial orbital semi major 
axis for the brown dwarf for the ef, = case. For each 
set of model parameters, we adopt four values for the 
brown dwarf's angle of apocenter, 0°, 90°, 180° and 
270° with respect to that of the binary system. 

In all our models, the orbital semi-major axes 
of the brown dwarf is larger than that of the bi- 
nary system. We are primarily seeking a critical 
initial semi- major axis (a c ), larger than which the 
brown dwarf survives the binary system's perturba- 
tion within a timescale Td- Our definition of survival 
is that the distance from the center of mass of the 
system to the brown dwarf (starting with all four 
values of the apoccntric arguments) must be smaller 
than a critical value Rd . The value of Rd is arbitrar- 
ily set to be 25 binary separations. In order to com- 
pare with the results of Holman & Wiegert (1999), 
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Fig. 1. The critical semi major axis as function of 
eccentricity when the mass ratio of the binary system 
fi — 0.1. The solid line represents the result for models 
with brown dwarfs been treated as massless test parti- 
cles which initially move on circular orbits. These models 
are analogous to those obtained by Holman & Wiegert 
(1999). The dotted line represents the results for mod- 
els in which brown dwarfs are assigned with finite mass 
fib = 0.05. The brown dwarfs are assumed to have circu- 
lar orbits about the center of mass of the binary system 
initially. The dash line represents the results of models 
in which brown dwarfs have mass [lb = 0.05 and have 
initial eccentricity e b = 0.4 

we choose Td = 10 4 binary period. Based on sev- 
eral test runs, we find that the value of a c does not 
change significantly if Td is increased to 10 6 binary 
period. Thus, we find a c to be a useful parameter 
to classify our results (Dvorak 2004). Although we 
use a totally different numerical scheme as Holman 
& Wiegert (1999), we are able to precisely reproduce 
the results in Table 7 of their paper when we set the 
mass of brown dwarf to be zero as they have. But in 
general, we choose [ib — 0.05. 

From these models, we find that the brown 
dwarf's "escape zone" (with semi-major axis a < a c ) 
is expanded slightly when they have finite mass (see 
Fig.l & 2 for the comparison.) The expansion of 
the "escape zone" is larger for /j, = 0.1 cases than 
that for the fi = 0.5 cases because the motion of 
secondary star is more affected by the finite mass 
of the brown dwarf. This effect is particularly no- 
ticeable for the fi = 0.1 and e* = 0.6 case where a c 
is expanded from 3.9 of massless particles to 6.0 of 
brown-dwarfs with /Xf, = 0.05. The hydrodynamical 
simulations indicate that fragmentation of circumbi- 
nary disks occurs primarily at around 1-2 binary sep- 



BROWN DWARFS 



3 



10 




Fig. 2. The same as in Fig.l except the mass ratio of the 
binary system /i = 0.5. 

aration away from their center of mass because this is 
the location where disk gas may accumulate as a con- 
sequence of binary star's tidal torque. Thus, most of 
the low-mass fragments formed in the circumbinary 
rings have a high probability of being ejected by the 
gravitational perturbation of their host binary sys- 
tems. 

We now consider a series of models with eb = 
0.4 while all other parameters are similar to those 
for the eb — case. The increases in a c in Fig. 1 
and Fig. 2 clearly indicate that brown dwarfs with 
eccentric orbits are definitely less stable than those 
with circular orbits. 

2.2. Large radial excursion of marginally stable 
systems. 

In general, the binary systems and the fragments 
formed in unstable circumbinary rings have non cir- 
cular orbits. Thus, most brown dwarfs formed close 
to the binary are likely to be ejected. But, brown 
dwarfs with initial semi major axis ab ~ a c , can be 
scattered to large distances from the center of mass 
of the system without escaping from its gravitational 
potential. 

We illustrate three such examples each with ab ~ 
a c . In model 1, we choose /i = 0.1, e* = 0.4, ab = 4.4, 
eb = and the argument of brown dwarf's apogee, 
9b = 90°. In this case, the brown dwarf reaches to 
200 binary separation by the end of simulation, i.e. 
t = T d . In model 2, (/z = 0.1, e* = 0.6, a b = 7.8, 
eb — 0.4, and 9b = 180°), the brown dwarf's orbit 
expands to 50 times the binary's initial separation 
at t ~ 0.6Td- But subsequently at t = Td, the extent 



of the brown dwarf's radial excursion is reduced to 
approximately its initial value. In model 3, (/i = 0.5, 
e* = 0.2, a b = 6.3, e b = 0.4, and 9 b = 270°), the ex- 
cursion reaches 100 initial binary separations. These 
examples indicate that wide and marginally stable 
orbits exist and that under some marginal circum- 
stances, brown dwarfs can be scattered to large dis- 
tance from but remain bound to some main-sequence 
binary stars. The recent discovery of a brown dwarf 
candidate at a distance of ~ 100 AU from a young 
binary star, TWA 6 Hyd may be examples of such a 
system. 

2.3. Ejection speed of escapers 

Brown dwarfs with ab < a c are ejected from the 
gravitational potential of their host binary system. 
We now examine their escape speed by series of cal- 
culations. 

From the distribution of the escape speed, we 
found that the ejection speed is typically half the 
orbital speed of the binary. For Rb ~ 3 x 10 15 cm, 
a ~ Rb, and the total mass of binary system ~ 1M Q , 
the binary's orbital speed would need to be ~ 3 — 5 
km s _1 . Our results show that the escaping speed 
of the brown dwarf ejecta is ~ 1 — 3 km s _1 . In a 
young stellar cluster, such as the Orion, this ejection 
speed is a fraction of the velocity dispersion of the 
cluster which is in a dynamical equilibrium. Thus, 
brown dwarfs ejected from the close proximity of the 
binary would not generally escape the gravitational 
potential of the cluster. This result is consistent with 
the large concentration of brown dwarfs in young 
stellar clusters such as the Orion complex (Lucas & 
Roche 2000). 

3. FORMATION OF BROWN DWARF PAIRS 

Indeed, close brown dwarf binaries have been 
found (Koerner et al. 1999), but these systems 
are generally not orbiting around some other binary 
main sequence stars. Similar to single brown dwarfs, 
close binary brown dwarfs may also be strongly per- 
turbed by the gravity of the binary and be ejected. 

3.1. Survival of pairs 

In order to test the survival probability of the 
brown dwarf binaries, we first place a massless test 
particle to simulate the dynamics of a secondary 
companion around the brown dwarf (For some in- 
teresting cases, a series of models with masses for 
both brown dwarf companions are included). This 
approximation allows us to first explore the range 
of parameters which may be favorable for the sur- 
vival of the brown dwarf pairs. Based on the results 
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in Fig.l and Fig. 2, we can identify the range of 
model parameters which leads to ejection of brown 
dwarf fragments. As a test, we adopt fi = 0.5 for 
the binary star and choose fib — 0.05, a& = 2.3 and 
apogee at 0° for the primary of the brown dwarf bi- 
nary. The secondary of the brown dwarf binary is 
assigned an initial semi major axis 2.5 and apogee 
at 0°. Thus the separation of the brown-dwarf pair 
is 0.2 which is inside the Roche radius of the primary 
(R R = (A*b/3) 1/3 = 0.25). We also assume that the 
center of mass of the brown-dwarf pair is initially on 
a circular orbit around the center of mass of the bi- 
nary system and the brown-dwarf secondary is on a 
circular orbit around the brown-dwarf primary. 

Then we try four cases of different eccentrici- 
ties of central binary stars relative to each other: 
e* = 0.0, 0.2, 0.4, 0.6. Our results indicate that 
the brown-dwarf pairs would survive their ejection 
from the neighborhood of the binary in the low- 
eccentricity (e* = 0.0 and e* = 0.2) limit. But, in 
the limit that the binary system has a large eccentric- 
ity (i.e. for the e* = 0.4 and e* = 0.6 models), the 
brown dwarf pairs have a tendency to become dis- 
sociated. Therefore it is plausible that brown dwarf 
pairs may remain bound to each other during their 
ejection from the binary system's gravitational po- 
tential. But it is also possible for their ejection to 
produce two freely floating single brown dwarfs. 

The above approximation provides a useful tool 
for us to identify the range of parameters which al- 
lows a brown dwarf pairs to survive. The mass- 
less approximation for the secondary is applicable to 
brown dwarf binaries with extreme mass ratios. We 
now takes the next iteration by replacing the (/it, = 
0.05) primary and a massless secondary brown-dwarf 
pair with a system of two equal-mass (fib = 0.025) 
brown dwarf companions. We find, with identical 
four sets of model parameters as above, brown dwarf 
binaries remain intact as they are ejected by their 
host binary stars. 

We also enlarge the separation of the brown- 
dwarf binary to 0.3, which is larger than its Roche 
radius. Again all four sets of initial conditions are 
used. In all these cases, the brown dwarf binary is 
always disrupted during its close encounters with the 
binary star. 

3.2. Pair capture 
We now explore the possibility that both brown 



dwarfs were formed as single objects and that they 
have captured each other to become a binary. In or- 
der to evaluate this probability, we repeat the earlier 
simulation in which 24 particles are placed in a ring 
around the binary system. The main difference from 
the earlier models is that a mass of fib = 0.05 is as- 
signed to each particles. The corresponding Roche 
radius for each individual particle is ~ 0.25 which is 
comparable to their initial separation. All the par- 
ticles are ejected from the gravitational potential of 
the binary system but no particle captured any other 
particle. We did not see any case in which the cap- 
ture happened. We thus conclude that this second 
scenario is very unlikely. 

4. SUMMARY 

For single brown dwarf satellites around binary 
systems, our dynamical calculations show that when 
they are formed in dynamically unstable regions, 
they are likely to be ejected from the gravitational 
potential of the binary system. These results pro- 
vide an explanation for the common sighting of field 
brown dwarfs. 

For binary brown dwarf satellite pairs, the 
calculations of four body (main sequence binary 
with brown dwarf pair) interaction show that these 
brown-dwarf pairs can remain to be bound to each 
other during the ejection if their initial separation 
is well within their Roche radius (which is typically 
Rr ~ 0.2 — 0.25 binary separation). 

We are grateful for friendly conversation with 
Matthew Bate, Michael Sterzik and Bo Rcipurth. 
Particularly, Michael Sterzik and Bo Rcipurth 
pointed out that the words "Main-Sequence" proba- 
bly should be replaced by "Solar- Type Stars" in the 
viewgraphs of the oral presentation of this paper. 
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